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SUMMARY 


Flight tests at Mach nvunbers from 0.8 to 1.4 were conducted on 
three rocket-propelled general research models of an airplane config- 
uration. Two of the models had thin unswept tapered wings of aspect 
ratio 3 and hexagonal airfoil sections. The two wings had different 
structural stiffness characteristics. The third model had no wing. 
Static and dynamic longitudinal stability, control, trim, and drag 
characteristics were obtained. The separate effects of the various 
airplane components on some of the characteristics were determined. 

The data were obtained by analyzing the response of the models to rapid 
deflections of the horizontal tail. 

The results obtained indicated some nonlinearity with lift coef- 
ficient of the lift-curve slope and static stability characteristics 
at high subsonic speeds. The variation of the lift-curve slopes with 
Mach number was fairly gradual. The variation with Mach number of the 
static stability was rather irregular, which was apparently caused by 
the addition of the wing to the conf ig\n:ation. The addition of the 
wing also caused the damping-in-pitch factor to vary considerably with 
Mach nimiber in the transonic region. Buffeting of the models with 
wings was observed at high subsonic speeds at lift coefficients a little 
below the maximum. The drag-rise Mach number was about O. 9 O for this 
configuration. The fuselage accoimts for about one-half of the mini- 
mum drag of the configuration. The variation of drag with lift coef- 
ficient indicated that the resultant aerodynamic forces on the wing 
and tall were approximately normal to the chord lines. 
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INTRODUCTION 


A general research program has been initiated by the NACA to 
determine by means of rocket -propelled models in free flighty the tran- 
sonic and supersonic longitudinal stability^ control and drag charac- 
teristics of airplane and missile configurations (references 1 to . 

The information is obtained by recording and analyzing the response of 
the models to intermittent disturbances in pitch applied to the models 
as they traverse the speed range. Presented herein are the results 
from the flights of two models of an airplane configuration having an 
unswept tapered wing of aspect ratio 3 QJ^d 4. 5-percent-thlck airfoil 
sections and of one similar model without a wing. From these results 
the total effects of the wing can be determined. The two models with 
wings were identical except that one had a solid-steel wing, whereas 
the other had a so lid -aluminum wing. The total effects on the config- 
uration of a change in wing flexibility thus appear in the results. 

Part of the data in a preliminary form and a detailed discussion of the 
analysis method have been presented in reference 1. 

For the models discussed in this paper, all-movable horizontal 
tails were used for elevator control, and during the flights the hori- 
zontal tails were deflected in alternate positive and negative direc- 
tions in approximate square-wave programs. The basic aerodynamic param- 
eters defining the longitudinal stability, control effectiveness, and 
drag characteristics were derived from the flight time histories for a 
Mach number range from 0.8 to 1.4. 

The models were flown at the Langley Pilotless Aircraft Research 
Station, Wallops Island, Va. 


SYMBOLS 


normal-force coefficient 



chord-force coefficient 



Cl lift coefficient cos a - sin a) 
Cp drag coefficient (C^ cos a + sin a) 
Cpjjiin minimum drag coefficient 


pitching -moment coefficient 
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normal acceleration as obtained from accelerometer, 
feet per second per second 

longitudinal acceleration as obtained from accelerometer, 
feet per second per second, positive forward 

acceleration of gravity, feet per second per second 

free-stream static pressure, pounds per square foot 

standard sea-level static pressure (2116 pounds per 
square foot) 

velocity, feet per second 


Mach nmber 

specific heat ratio (l. 40 ) 

Reynolds number, based on wing mean aerodynamic chord 

wing area (including the area enclosed within the fuselage), 
square feet 

horizontal tail area (including the area enclosed within the 
vertical tail) , square feet 

lateral distance from side of fuselage, inches 

wing semispan measured from side of fuselage, inches 

wing mean aerodynamic chord, feet 

tip chord of wing, inches 

torsional modulus of elasticity, pounds per square inch 

wing torsional-stiffness parameter, inch-pounds per radian 

couple applied near wing tip in plane parallel to model 
center line and normal to chord plane, inch-pomds 

local wing twist angle produced by m measured in plane 

parallel to that of m, radians (when used as a subscript) 


dynamic pressure 



4 


NACA RM L 5 OKI 6 


e 


a 

S 


€ 



angle of pitch, degrees 

moment of inertia about y-axis, slug-feet^ 
angle of attack, degrees 
elevator deflection, degrees 
dovnvash angle, degrees 

dovnvash angle at zero angle of attack, degrees 
period of pitching oscillation, seconds 
time, seconds 

time to damp to one -half amplitude, seconds 

pitching -moment coefficient at zero angle of attack and zero 
elevator deflection 


Subscripts: 

T trim 


a 


_ da c 

~ dt W 


q = ^ -L 

dt 2V 


W ving 


f fuselage 

The symbols a, a, 5, and q used as subscripts indicate 
the derivative of the quantity vith respect to the subscript, for 



example^ 
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M0DEI5 AND APPARATUS 
Models 


A drawing of the complete airplane configuration investigated is 
shown in figirre 1. The wingless model is shown in figure 2. For the 
wingless model, it was necessary to add a vertical tail to the bottom 
of the fuselage to prevent large rolling motions resulting from yaw 
disturbances. Photographs of the models are shown in figure 3 . 

The model was designed as a general research vehicle for investi- 
gating the effects on airplane stability, control effectiveness, and 
drag characteristics of changes in the various airplane components. 

The fuselage was not an optimum shape for low drag but was selected 
from considerations of ease of fabrication and adaptability for 
altering various components. The fuselage was composed of a cylindrical 
center portion with nose and tail sections defined by the ordinates in 
table I. A fairly large vertical tail was provided to insure adequate 
directional stability for widely different wing configurations. The 
following letter symbols are used throughout this paper to designate 
the three models: 


Model 

designation 

Description 

A 

Steel ving 

B 

Aluininuin wing 

C 

Wingless 


The wings and horizontal tails had 4. 5-percent-thick hexagonal 
airfoil sections as shown in figure 1. Model A had a solid- steel wing 
(reference l) and model B had a solid-alumin\im wing. The horizontal 
tails were solid al\amimm on all three models. 

The horizontal tail was mounted on a ball bearing built into the 
vertical tail and was rotated about a hinge line located at 42 percent 
of the tall mean aerodynamic chord (fig. 4). The horizontal tail was 
deflected in alternate positive and negative directions during the 
flights by means of an electric motor and cam arrangement in models A 
and B and a hydraulic control system in model C. The control-position 
pickup was located on the lower end of the elevator push rod. The 
measvired elevator motion can therefore be considered as the motion of 
the root of the elevator. As can be seen in figure 4, a gap existed 
between the vertical tail and the root of the horizontal tail. This 
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gap extended over the foivard and aft 30 percent of the horizontal tail 
root chord. Its vidth was approximately l/2 percent of the chord. 

An additional noninstrumented model shown in figure 5 was flown to 
determine the drag of the extra vertical tail used on model C. Two 
scaled-down vertical tails were moimted as wings on a f in-stahilized 
body of revolution and flown at zero lift to determine the minim\mi drag. 
The model drag was determined from the velocity-tim.e curve obtained by 
the Doppler radar. The drag of the vertical. tails was obtained by sub- 
tracting the drag of the body plus fins (reference 5 ) from the drag of 
the complete model. 


Instrumentation 

Models A and B contained six-channel telemeters transmitting 
measurements of noimial and longitudinal acceleration, elevator deflec- 
tion, angle of attack, total pressure, and a reference static pressure. 
Model C had a seven-channel telemeter giving records of the same six 
quantities plus an additional static pressure. Angle of attack was 
measured by a vane -type instrument located on the nose of the model 
(figs. 1 to 3)- This device is more fully described in reference 6. 

The total-pressure tube was located on a small strut below the fuselage 
(figs. 1 to 3)* On all three models the reference static pressure 
measured was the pressure inside the cone of the angle- of- attack vane. 
Model C also had a static-pressure orifice located on the top of the 
fuselage at 4.9 inches behind the beginning of the cylindrical section. 
The total- and static-pressure locations had been calibrated previously 
for zero angle of attack on instrumentation test models. 

Doppler radar and tracking radar units were used for obtaining 
checks on model velocity, range, and elevation as functions of time. 
Atmospheric conditions were determined from radiosondes released shortly 
after each flight. Fixed and manually operated l6-millimeter motion- 
picture cameras were used to photograph the launching and the first 
part of the flights. 


Launching Procedure 

The models were boosted to maximum velocity by a 6 - inch -diameter 
solid-fuel Deacon rocket. The models contained no sustaining rockets. 
The tapered after ends of the models fitted inta adaptors on the booster 
rockets (fig. 6 ) . This type of attachment permitted freedom in roll and 
fo27ward movement of the model with respect to the booster but restrained 
the model in yaw and pitch. Separation of the models from the boosters 
was accomplished by reason of the different drag/weight ratios of the 
m.odel and booster following booster burnout. The models were laimched 
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at angles of approximately 45^ from the horizontal, by means of a 
crutch-type launcher as shown in figure 6. 


The measured mass characteristics of the models are given in 
table II. The torsional stiffness characteristics of the wings were 
determined by applying a couple at the tip and measuring the twist at 
several sections. The results are shown in figure 7* The factor 


independent of size or material of construction. Measured points from 
three other geometrically similar wings, one larger (unpublished) and 
two smaller (reference j) than those described herein, are included 
on the curve in figure 7. For use in comparing the aeroelastic prop- 
erties of these wings with other results, the values of free-stream 
static pressure divided by the standard sea-level pressure are shown 
in figure 8 for each of the three models. 


The data on the characteristics of the models were obtained during 
the decelerating part of the flights following separation of the model 
and booster. The Doppler radar obtained velocity data for all three 
models dtiring the boost period but failed to track the models after 
booster separation. The tracking radar obtained flight-path informa- 
tion for all models during the boost period and for model C during 
the entire flight. For models A and B the Mach numbers and dynamic 
pressures during decelerating flight were therefore determined from 
the telemetered total and static pressures. 

The angles of attack measured by the vane on the nose of the model 
were corrected to angles at the model center of gravity by the method 
of reference 6. 

The Reynolds numbers (based on the wing mean aerodynamic chord) 
obtained during the flights are shown in figure ^ as a function of Mach 
number . 

The models discussed herein served as development models for 
evaluating and improving the testing and analysis procedures, and as 
such, incorporated various changes as the program proceeded. For 
model A, elevator deflections of approximately ±2^ were used with a 


Preflight Measurements 


plotted. 



is a nondimens ional parajneter which makes the result 


TEST AND ANALYSIS PROCEDURES 


Tests 
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programming rate of one cycle per second in an approximate square-vave 
pattern. This rate vas too fast to permit the determination of damping 
or trim characteristics at the lover Mach numbers (reference l) . The 
elevator on model B was deflected between + 4 ^ and - 2 ° at a rate of about 
one cycle in I.7 seconds and the center of gravity was placed farther 
forward to give smaller oscillation periods. 

During the first elevator deflection (at M = 1 . 37 ) following 
booster separation on model B the elevator moved about 1 ^ and then 
stopped^ probably because the hinge moments were too great for the power 
system. At a Mach number of about I.I7 the elevator began moving again 
and operated normally thereafter. Between these Mach numbers no oscill- 
ation data were obtained. At subsonic speeds the elevator did not remain 
at a fixed deflection for either model A or B during some of the oscil- 
lations. The elevator was overbalanced at subsonic speeds and apparently 
moved off the stop due to accidental play in the system as the hinge 
moment reversed during the oscillation. For model A the amount of 
elevator motion was small (0.25^ or less) and its effect on the model 
motion was negligible. For model B the elevator motion was about O.60 
and caused an appreciable change in the character of the oscillation. 
These oscillations were not used in the determination of stability and 
damping data. 

Because of the difficulties encountered on models A and B a new 
control system was designed for model C. This was a hydraulic system 
having a much greater power output than the original electric motor and 
cam system^ and having a positive stop for the elevator. The deflection 
of the elevator was more nearly a step function than for the other two 
models and the elevator remained at a fixed deflection during all of the 
free oscillations. The elevator on model C was deflected between +30 
and -1.2° at a rate of one cycle in I.9 seconds. 


Analysis 

After each elevator deflection the models experienced short-period 
transient oscillations of angle of attack and normal and longitudinal 
accelerations. These oscillations and the subsequent steady-state values 
were analyzed by the method described in appendix A of reference 1 to 
find the longitudinal aerodynamic characteristics. The method is essen- 
tially a graphical procedure with a minimum of computational work. Its 
success depends to a large extent upon the accurate measurement of 
angle -of -attack changes of the oscillating aircraft to permit a direct 
determination of the lift-curve slope and the control effectiveness. 

The method does not require the reduction of the data to frequency- 
response form before the determination of the aerodynamic derivatives. 

If desired^ the frequency-response characteristics can of course be cal- 
culated from the aerodynamic derivatives. 
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All of the procediires so far proposed for determining aerodynamic 
derivatives from free -flight tests assume that the aerodynamic coeffi- 
cients are linear functions of the variables. The results shown herein 
and other unpublished results of rocket model tests indicate that in the 
transonic region this will often not be true. The nonlinearities usually 
appear as different slopes Ct and C measured when the model is 

oscillating over different ranges of lift coefficient. When nonlineari- 
ties appear, the slopes obtained should be considered as average values 
over the range of lift coefficient covered. Reference 8 contains a 
detailed study of the effects of nonlinear derivatives on aircraft 
motions. Considerable Judgment is sometimes necessary in interpreting 
the results obtained when nonlinearities exist. For example, reference 1 
lists three ways by which can be found. If the aerodynamic 


quantities are nonlinear, none of the methods give a true answer. 
However, experience with a number of rocket models shows that the method 
of determining Cl^ which is most nearly correct depends on the config- 

ation tested. For an airpleine configuration with a fairly conventional 
horizontal tail the distance to the center of pressure of the lift due 
to control deflection can be estimated fairly accurately and the most 
accurate value of will probably be that obtained by dividing 

^mg this distance expressed in terms of the mean aerodynamic chord. 


For a missile configuration having all-movable wings near the center of 
gravity for control, this distance cannot be accurately estimated. For 
such a configuration C^g will be of the same order of magnitude as 

C]^ and thus equation (A-I7) in reference 1 relating these quantities 


to Clr should give a more nearly correct result. 

°trim 


For determining the period and damping of the oscillations , the 
angle-of -attack time history is used rather than the lift coefficient. 

The equations show that the angle-of-attack oscillation is affected only 
a small amount by the lift- curve slope, whereas the lift coefficient 
oscillation is directly dependent on the lift-curve slope. Thus, if a 
nonlinear C-^ exists, the period and damping obtained from the angle- 

of-attack oscillation are more nearly a correct measure of the static 
stability and damping characteristics (see reference 8). The errors 
involved in certain simplifying assumptions made in the analysis are 
examined in some detail in appendixes A and B of reference 1. 

The models discussed herein were not roll- stabilized and were 
thus free to roll under the action of any structural asymmetry. The 
actual rate of roll during flight was not measiired. Reference 9 shows 
that the steady rate of roll required to alter the longitudinal stabil- 
ity significantly is of the same order of magnitude as the natural 
frequency in pitch. For the lowest pitching frequency encountered on 
these models this would necessitate a rate of roll of about 3 revolutions 
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per second. The differential wing incidence required to produce this 
motion would be about 2 ^ or greater. Since the wings were machined in 
one piece^ the actual differential incidence should be less than 0 . 1 ^. 
The effects of roll on the data presented herein should therefore be 
negligible. 


ACCURACY 


It is impossible to establish the absolute accuracy of the measured 
quantities on any given model because the instrument calibrations can 
not be checked during or after the flight. The repeatability of data 
obtained from identical models or from identical components on different 
models furnishes the best check on the over-all accuracy. In general, 
the possible instrument errors should be proportional to a certain 
percentage of the total calibrated range of the instrument. Most of 
the probable instrumentation errors occur as errors in absolute magnitude. 
Incremental values or slopes and the variation with Mach number of the 
various quantities should in general be more accurate than the absolute 
values. The following table gives estimated values of the possible 
systematic errors in the absolute values of and as eiffected 

by the accelerometer calibration ranges: 


Model 

M 

ACl 

ACd 

A 

0.80 

±o.o 4 o 

± 0.020 

B 

.80 

i.Ol+O 

±.005 

C 

.80 

±. 0 l 4 

±.005 

A 

1.30 

±.015 

±.008 

B 

1.30 

±.015 

±.002 

C 

1.30 

±.005 

±.002 


The random errors may be judged by the scatter of the data points 
shown in the figures. Model A was the first model flown in this 
research program and the longitudinal accelerometer was calibrated to 
cover the thrust range as well as the drag range, to study the perform- 
ance during the boosted phase of flight. The resulting accuracy for 
measuring drag was not considered satisfactory (see table above) so the 
longitudinal accelerometers were calibrated to cover only the drag range 
on models B and C. 

Further errors in the aerodynamic coefficients arise because of 
possible dynamic-pressure inaccuracies which are approximately twice as 
great as the errors in Mach niomber. A consideration of all the factors 
involved indicates that the Mach numbers are probably accurate to about 
±1 percent at peak velocity when Doppler radar information is obtained. 
The Mach numbers for all models are probably accurate to ^ percent or 
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better at supersonic speeds. At subsonic speeds the Mach numbers are 
probably accurate to ±2 percent on models for which flight-path data 
are obtained by radar (model C) and are less accurate for models on 
which radar tracking is not obtained (models A and B) . A very valuable 
check on the accuracy of the supersonic Mach numbers has been found to 
exist in the static^pressure measurements. The static -pres sure sources 
that have been calibrated on rocket models show an abrupt pressure change 
at or near a Mach number of 1 . 0 . If the Mach number at which this 
change occurs has been accurately determined on previous models by means 
of Doppler and flight-path radar, then the Mach number at this point on 
subsequent flights is known irrespective of errors in the pressure 
measurements. By using the Doppler radar at peak velocity and this 
pressure change at a known Mach number, the pressure measurements may be 
corrected for any zero shifts that may have occurred. 

The errors in the measured angles of attack and elevator deflections 
should not vary with Mach number because they are not dependent on 
dynamic pressure. Probably the greatest possible error in angle of 
attack is caused by possible aerodynamic asymmetry of the angle-of- 
attack vane which is not detectable prior to flight. The elevator 
deflections should be accurate to about ±0.1^ and the increments in angle 
of attack to about ±0.2^. 


RESULTS ARD DISCUSSION 


As noted in table II, the models were flown with different center- 
of -gravity locations, but the aerodynamic derivatives for model B 
presented in the figures have all been converted to values for a 
center-of -gravity location of 12.^ percent of the mean aerodynamic 
chord for comparison with models A and C. Because of the nonlinearity 
with lift coefficient of some of the aerodynamic derivatives at sub- 
sonic speeds, the conversion of the data to a different center of 
gravity may lead to some decrease in the accuracy of the results. 


Lift 

Figure 10 shows typical plots of the lift cuirves obtained for the 
three models. The data points shown are those reduced from the telem- 
eter records at time intervals of 0.02 second. The data were recorded 
continuously and any number of points could be obtained for such a plot. 
The data obtained were plotted as lift coefficient against angle of 
attack during the first 1^ or 2 cycles of each oscillation. 

The data for model C in figure 10 are plotted to the same scale as 
the data for models A and B. As pointed out in the section "ACCURACY,” 
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however, the lift coefficients for model C are considerably more accurate 
than for the other two models because of the more sensitive accelerometer. 
In the process of finding lift-curve slopes the ciorves for model C were 
plotted to a much larger scale commensurate with their accxiracy. Two 
ciorves are faired through each set of data illustrated for model C in 
figure 10, 'one curve connecting points measured while the angle af attack 
was increasing and one curve for points obtained while the angle of 
attack was decreasing. The difference is in such a direction as to appear 
as a phase lead of the lift coefficient ahead of the angle of attack, 
which is in the opposite direction to what mi^t be expected from poten- 
tial-flow unsteady -lift considerations (reference 10) . The same effect, 
of approximately the same magnitude, can be detected in some of the 
curves for models A and B when plotted to the larger scale used in the 
data analysis. In some cases the scatter of the data points for 
models A and B, attributable to the less sensitive accelerometers, masks 
the difference cavised by the direction of angle- of- attack change. On 
models A and B the effect is much larger at subsonic speeds than at 
supersonic speeds and is larger at high lift coefficients than at low lift 
coefficients. 

Part of the difference in lift at a given angle of attack can be 

accounted for by the lift derivatives C^. and Ct . The calculated 

a q 

effect from estimated values of Ct . and Ct is shown for several 

■^a q 

curves in figure 10. The largest difference in lift at a given angle 
of attack occurs for model A at the points where the airfoil appears to 
have stalled. The sajne effect on stalled airfoils is observed in static 
wind-tunnel testing. The lift data in reference 4 show the same effect 
of rate of angle -of -attack change on the lift, but the effect is greater 
in reference 4 than for the models described herein. The model in 
reference 4 also had a wing with hexagonal airfoil sections, the rear 
wedge of the airfoil being more blunt than those on the models herein. 
Reference presenting data for a model with a circular-arc airfoil, 
showed only a very small effect of rate of change of angle of attack. 

The possibility exists that the effect under discussion is caused by 
asymmetric air-flow separation even at low lift coefficients and is a 
function of the type of airfoil section. 

The lift-curve slopes obtained for the three models are shown in 
figure 11. Some additional points from small-amplitude oscillations are 
shown for model A which were not' included in the preliminary analysis 
of reference 1. Slopes measured from motions of small amplitude could 
be expected to be less accurate than those from larger amplitude motions. 

No large or abrupt changes in lift-curve slope with Mach n-umber occur 
in the transonic region. This behavior is probably attributable to the 
low aspect ratio and the thin airfoil sections of the wing and horizontal 
tail. Thicker airfoil sections or higher aspect ratios have been shown 
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to cause larger variation in lift-curve slope at high subsonic speeds 
(reference ll) . 

The data from all three models show some nonlinearity of the lift 
curves with angle of attack in the transonic region. The relative magni- 
tude of the nonlinearity is greater for the wingless models probably 
because of the greater proportionate effect of the fuselage. 

Differences occurred between the slopes measured for the two winged 
models (fig. 11(a)), model B having higher slopes up to M = I.I 5 and 
lower slopes at M = l.U. A calculation of the effect of wing twist at a 
Mach number of 0.8, using the stiffness curve in figure 7^ showed a 
possible difference in of 2 percent which is within the accuracy of 

measuring the slopes. The increase in dynamic pressure with Mach number 
and a movement of the wing aerodynamic center from a position forward of 
the elastic axis to one aft would account for the differences in lift- 
curve slope shown. 

Figure 12 shows the lift-curve slopes of the tail and of th^ wing 
plus interference near zero lift. The values for wing plus interference 
were obtained by subtracting the slopes for model C from those for 
model A. The rocket model data for the wing lift-curve slope also 
include a negative increment in slope caused by the wing downwash acting 
on the tail. The wing-alone slopes would therefore be somewhat higher 
than those shown in figure 12. Also shown are values of lift-curve slope 
of the wing alone from unpublished tests on the transonic bump in the 
Southern California Cooperative Wind Tunnel of a similar wing and 
theoretical results at supersonic speeds for a straight rectangular wing 
having an aspect ratio of 3*0 (reference 12). The wind-tunnel and 
rocket model data show similar variations with Mach number. The agree- 
ment between the rocket model results and the theoretical calculations is 
relatively good. Since the lift added by the wing on the rocket model 
is as great or greater than that for the wing alone as tested in the 
wind tunnel, the indications are that the wing area Included in the 
fuselage (about 27*5 percent) may be considered almost fully effective 
in producing lift on this configuration. Since the wing and horizontal 
tail are identical in plan form and airfoil section, the lift effective- 
ness of the tail (fig. I 3 ) should equal the lift-curve slope of the wing 
alone when multiplied by the area ratio of the two surfaces. Such a 


the tail are probably caused by the gaps at the root of the tail. 

The values of Z^Cl/AB in figure I 3 are in good agreement for the 
three models. Because the tail lift was a much greater proportion of 
the total lift for model C, the values of ACl/AB are more accurate for 
this model than for the other two models. The tail lift-ciirve slopes 
should also be nonlinear with lift coefficient but the method of 



shown in figure 12. The smaller slopes shown for 
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determining /SCjjAb does not permit the determination of this effect 
and the data are thus shown as incremental values rather than slopes. 

If it is assumed that the dc/da acting on the tail of the wing- 
less model is negligible, then the difference between the lift-curve 
slope Cx^ for model C and the ACl/A 6 for model C will give the 
lift-curve slope of the fuselage. The nonlinearity of the curves of 

prevents this procedure from being used, but at the lowest and 

highest Mach numbers attained the indications are that the fuselage 
slope is of the order of 0.002 or O.OO3, which approaches the accuracy 
of measuring the slopes on this model. 


Maximum Lift and Buffeting 

The data on models A and B indicated that maximum lift coefficients 
and some buffeting information were obtained over a range of high 
subsonic speeds. A discussion of the maximum lift coefficients is given 
in reference 1 and the buffeting information is discussed in somp detail 
in reference I3. Figure l4 is a summary of the information obtained 
showing the maximum lift coefficients and the region where buffeting 
was observed. A curve is also shown defining the highest lift coeffi- 
cients reached. At these values the maximum lift coefficient of the 
configuration has not yet been reached. Below this cujcve and outside 
the cross-hatched area no evidence of buffeting appeared in the records. 


Drag 

Minimum drag coefficients for all three models are shown in 
figure 15. The data for model A, from reference 1, have been omitted 
below a Mach number of I.07 because they are inaccurate at subsonic 
speeds as indicated by a comparison of the measured minimum drag coeffi- 
cients for model A with those for models B and C, which had considerably 
more sensitive longitudinal accelerometers. At Mach numbers above 1.0 
the accuracy is much better and the agreement between models A and B 
was good. 

The drag divergence Mach number for this configuration is apparently 
about 0.90. A breakdown of the minimum drag coefficient into its 
various components is shown in figure I6, all coefficients being based 
on total wing area. The drag of the vertical tail was determined from 
the special model flown for this purpose. The drag for fuselage-plus - 
tail surfaces was obtained by subtracting the drag of one vertical tail 
from the drag of model C. This drag subtracted from that for models A 
and B gives the drag of the wing. Multiplying the minimum drag of the 
wing by the ratio of exposed horizontal tail area to exposed wing area 
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gives the drag of the horizontal tail. The drag of the fuselage vas 
found by subtracting the drag of two vertical tails and one horizontal 
tail from the drag of model C. The drag values obtained as described 
of course include some unknown interference effects. The fuselage drag 
is seen to be about one-half of the minimum drag of the airplane at all 
Mach mmibers. When this drag is compared with that of the parabolic 
bodies tested in reference l4^ it is evident that the fuselage drag 
is fairly high when based on maximum cross-sectional area but is not 
unreasonable when based on the fuselage volume. 

The minimum drag coefficients of the wing and horizontal tail can 
not be considered as very accurate when obtained in the manner described. 
For comparison, the wing drag, obtained from the data of reference 7 aj^d 
unpublished body drag data, is shown in figure I 6 . The wings on the 
models in reference 7 were similar to those on models A and B, but they 
had deflected ailerons and were rolling continuously. The wing minimum 
drag coefficient, obtained as described herein, appears to be of the 
right order of magnitude and is only a small part of the minimum drag of 
the configuration. 

The variation of drag with lift is shown in figure 17 in the form 
of dCpy^CL^. Although the minimum drag for model A was inaccurate at 
subsonic speeds, it appears that the factor is of the right 

order of magnitude when compared with model B although the scatter of 
the data points is greater than for model B. This is in accord with 
the statements in the section on "ACCURACY" that most of the errors in 
telemetered quantities occur as drifts in the zero value rather than 

changes in slope. Comparison of the values of dCjy^C^^ shown for 
models A and B with the factor i/57.3CL(x i*or both models indicates 
that the resultant force is approximately normal to the wing chord. 
Comparison of the values of dCp^dCL^ for model A at M = 1.35 with 
those for model B indicates that the values for model A are somewhat 
larger. The reason for this difference is not definitely known but the 
values for model A were obtained from oscillations of larger amplitude 
and the possibility of some nonlinearity of dCiy^dCL^ with lift 
coefficient exists. 

Values of meiximum lift-drag ratio and of C^ at which maximum 

lift-drag ratio occurs are shown in figures I 8 and I 9 . The symbols 
designate points where the maximum values were actually measured. At 
other Mach numbers the values represented by the faired curve were 
obtained by an extrapolation of the drag curves based on the measured 

values of ^Dmin dCj^/dCj^ . The relatively low values of (L/D)max 
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and high values of Cl for (L/D)max both a reflection of the high 

drag of the fuselage. 


Static Stability 

The measured periods of oscillation of the angle of attack are 
shown in figure 20. Not many points are available for model B because 
of the difficulties encoimtered with the control system as mentioned 
previously. The data converted to for a center of gravity at 

12.4 percent of the mean aerodynamic chord are shown in figure 21. 

Above a Mach mmiber of 1.0 some difference existed between models A 
and with model B exhibiting the greater stability. Further comments 
on this point will be made subsequently in the discussion of the trim 
characteristics. Model A indicated that Cj^ for the complete airplane 

varies with lift coefficient below M = O.92. The wingless model showed 
nonlinearity at Mach numbers above 0.86. 


Both models A and B showed the same type of variation of Cm^ with 

Mach number in the region where data were obtained for both models. The 
somewhat irregular variation with Mach number of the Cm^ for models A 

and B is apparently an effect of the wing on the configuration because 
Cin^ for model C varies smoothly with Mach n\mber. Subtracting the 

values of for the winged and wingless models gives the total 


effect of the wing on the stability of the airplane. The result is 
shown in figure 22 for the more rigid wing. For the center- of- gravity 


position used the contribution of the wing itself to ^ 




should be 


negative and the increment of AC 




caused by the wing downwash on the 


tail will be positive. It is not possible to tell from the data what 
the relative magnitudes of the two effects are. It is estimated that 
the ACiji^^ due to the wing alone would be about -0.01 at subsonic 

speeds and about -0.02 at supersonic speeds. 


By subtracting the of the horizontal tail from the Cjj^^ 

of the wingless model the stability contribution of the fuselage ACm^^^ 

is obtained. The result is shown in figure 22. This procedvire assvunes 
that the downwash over the tail due to the fuselage does not change with 
angle of attack. Because of the high tail position, this is probably a 
fairly good approximation. The value of AC^^^ calculated from 

reference 15 is O.OI7 which shows good agreement with the measured valVies. 
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The aerodynamic -center locations for the winged models A and B 
are shown in figure 23 . Both models show the same trends with Mach 
number in the region where data were obtained for both models. The 
aerodynamic center for model B was farther rearward than for model A 
at Mach numbers above 1.0. The large values of dynamic pressure occurring 
during the tests probably exaggerate the effects of wing flexibility 
compared to the probable operating conditions for a full-size airplane 
at the same Mach number. 

The aerodynamic center for the wingless configuration varies around 
an average position at about 200 percent of the mean aerodynamic chord 
or about 70 percent of the fuselage length from the nose. The fuselage 
instability causes an average forward shift in aerodynamic center of 
about 18 percent of the mean aerodynamic chord on the complete 
configuration. 


Damping in Pitch 

The time required for the oscillations to damp to one -half amplitude 
is shown in figure 24(a) and the data converted to the damping factor 
Cmq + and corrected to a center of gravity at 12.4 percent of the 

mean aerodynamic chord for model B are shown in figure 24(b) . The 
wingless configuration shows a damping factor that is nearly constant 
with Mach number and exhibits some nonlinearity with lift coefficient 
at subsonic speed. The winged configurations have damping factors that 
are of roughly the same order of magnitude as those for the wingless 
configuration but show considerable variation with Mach number, 
particularly between M = O. 9 O and M = 1.05* Although the damping 
factors Cniq + Cjn^ are approximately equal for the winged and wingless 

models, the actual damping of the oscillations is considerably better for 
the winged models (fig. 24(a)) because of the influence of the lift- 
curve slope on the damping. This effect would occur on full-scale 
airplanes, the relative effect of the lift-curve slope depending on the 
radius of gyration ( see reference l) . 

Theoretical calculations of the damping of rectangular wings at 
supersonic speeds are available in reference 12. The damping derivatives 
Cffig ^m^ were calculated for the wing alone and tail alone at 

Mach numbers of 1.25 and 1.40 by use of this reference. For the wing 
and tail combination on models A and B the increments in the derivatives 
due to the effect of the wing on the tail were added to the wing-alone 
and tail-alone derivatives. A large negative increment in 0^^ occurs 

because of the downwash lag effect and a small positive increment in 
Cm exists because of the downwash angle at the tail caused by the 
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wing lift derivative both calculations a value of 

was assumed. The calculations indicated that at M = 1.25 


^ = 0.40 
da 

the 


of the wing alone was positive and larger than the negative Cm^ of the 

wing, resulting in a dynamically unstable wing at M = 1.25 for the 

center of gravity of the tests. When the tail was placed behind the 

wing the negative caused by the downwash lag overcompensated for 

the positive contributions of the tail alone and wing alone to 


resulting in a negative for the complete configuration. 


The calculations also showed that the of the tail alone was 

small compared to the of the tail. The measured damping factor 

^m ^m* fo^ model C is therefore essentially the of the tail, 

q a q 

which is the largest contribution to the pitch damping factor and which, 

according to figure 24(b) , does not vary much with Mach number in the 

range covered by the tests. Since the measured lift-curve slopes of the 

tail are less than the theoretical, it might be expected that the measured 

damping would be less than the theoretical. Figure 24(b) indicates this 

to be the case for models A and B, but not for model C. The measured 

values for all the models of course include the damping contribution of 

the fuselage which is not included in the theoretical calculations. 


The work in reference 12 indicates that the positive values of 
for fhe wing alone and tail alone increase with decreasing Mach 

number at supersonic speeds. Also, the variations in static stability 
in the transonic region, figure 21, indicate that large changes in 
downwash are probably occurring. Thus it is probable that the variations 
with Mach number of the factor for the winged configurations 

are caused by the term which could become positive under conditions 

of small downwash. 


The differences in between models A and B are not 

necessarily attributable to wing flexibility but are more probably a 
measure of the experimental accuracy of determining the damping derivatives. 
The primary reason is that for rocket models similar to those under 
consideration the contribution of the Ct term to the damping is very 


large and the factor is the difference between two large 

quantities. Small errors in or in ^±^2 therefore have a 

large effect on the values of For models A and B the Cl^ 

term accounted for about two-thirds of the total damping, while for 
model C it accounted for about one-third. Thus, the damping derivatives 
for model C are probably numerically more accurate than for models A and B. 


NACA RM L5OKI6 


19 


Longitudinal Trim and Control Effectiveness 


The trim lift coefficients and angles of attack are shown in 
figure 25. The trim lift coefficients for model C are not shown because 
they are very small (about ^ 0 . 03 ) compared to those for models A and B 
and have no great significance. The complete configuration (models A 
and B) does not appear to have any large trim changes in the transonic 
region. Model C indicated a more or less abrupt trim change in both 
lift coefficient and angle of attack between M = O.9O and M = 1 . 0 . 

This same type of change may or may not have occurred for the complete 
configuration; the number of trim points obtained in this region was 
not sufficient to establish definitely the correct fairing of the curves. 


It would be expected that the curves of and oor for 6 = 2.0^ 

would coincide for models A and B in the Mach number region where data 
were obtained for both models. The curve for & = 2 . 0 ^ for model B 
was obtained by a linear interpolation from the measured values for 
b = - 2 . 0 ^ and 4 . 0 ^. The failure of the curves for models A and B to 
coincide may be due to an error in measuring the elevator deflection 
for either or both of the models or may be an indication of some non- 
linearity in the variation of ojj. and Ct with elevator deflection.* 


Both the cxip and curves for model B indicate smaller 


increments due to elevator deflection and ZsC 




than for 


model A in the region near a Mach number of 1.1 where the trim data 
for the two models overlap. Since the Ao^^AS is equal to 

and the tail surfaces were identical on the two models, the disagreement 
between the two models is an indication of greater stability on model B. 
This agrees with the aerodynamic -center data in figure 23. Since the 
two indications of stability were obtained independently of each other 
in the reduction of the data, it appears that the greater stability on 
model B at M = 1. 1 actually existed and is not the result of errors 
incurred in the data analysis. 


The pitching-moment effectiveness of the elevator is shown in 
figure 26. The all-movable tail provided an effectiveness which is 
practically constant with Mach number. It is probable that the dCm/dB 
of the elevator is nonlinear with elevator lift coefficient or deflection, 
but this could not be determined from the results because only two 
elevator deflections were used on each model. The values shown in 
figure 26 are thus shown as incremental values ACjy^lb and represent 
average values over the elevator deflection range covered by the tests. 

For this reason, agreement of the data for the three models is considered 
good. 
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Figure 27 presents the effectiveness of the elevator in changing 
the trim lift coefficient of the airplane. Since ACnyAS remains 
fairly constant with Mach number the decrease in with 

increasing Mach nmber is caused by the increase in static stability. 

The difference between the curves for the two models near M = 1.1 
again indicates the greater stability of model B, 

The pitching-moment coefficients at zero angle of attack and zero 
elevator deflection are shown in figure 28. The magnitude of is 

about the same for the winged and wingless models, which indicates that 
■the Cjn^ is primarily due to the fuselage-tail combination. The 

converging rear portion of the fiiselage (fig. l) would cause a downflow 
over the tail and the drag of the horizontal and vertical tail siirfaces 
would cause a positive pitching moment. From the drag data in figvire 16 
this latter effect is calculated to be = 0.010 at M = I.15 and 

Cirio = 0.006 at M = 0.85, which are small compared to the values shown 

in figiire 28. Since the wing and fuselage are symmetrical about the 
longitudinal axis the lift and thus the downwash due to lift should be 
zero at zero angle of attack. The elevator deflection required to trim the 
model at zero angle of attack should therefore be equal to the downwash 
at the tail which is producing the major part of the These values 

of downwash angle are shown in figure 29. Apparently this downwash angle 
is caused almost entirely by the fuselage, the angle being approximately 
1*^ at subsonic speeds and 2° at supersonic speeds. The different 
variations with Mach number of the Cuiq and Cq curves for models B 

and C in the region between M = O.9O and M = 1.00 is not necessarily 
attributable to the wing. As stated previously, the trim data for model B 
(from which Cjj,^ and eo determined) were not complete enough in 

this Mach number region to establish the correct fairing of the curves. 
Since the and €o are apparently caused mainly by the fuselage 

shape, it is probable that a variation with Mach number such as that shown 
for model C between M = O.9O and M = 1.00 would be obtained for the 
complete configuration if more trim data were available. 


CONCLUSIONS 


Flight tests of rocket-propelled models (with and without wings) 
of an airplane configuration having thin unswept tapered wings and 
horizontal tails of aspect ratio 3 indicated the following conclusions: 
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1 . The variation of lift coefficient vith angle of attack vas 
nonlinear in the transonic region. The variation of the lift-curve 
slopes with Mach number was fairly small and gradual. 

2. The portion of the wing area enclosed within the fuselage may 
be considered fully effective in producing lift on this configuration. 

3. An apparent reduction in lift effectiveness of the all-movable 
horizontal tail as compared to that of the wing occurred throughout 
the speed range; this reduction was probably caused by chordwise gaps 
at the root. 

Buffeting of the model occurred at high subsonic speeds at 
lift coefficients a little below the maximum lift coefficient. 

5. The minim-urn drag of the fuselage was about half of the drag 
of this configuration. The drag-rise Mach mmiber was about O.9O for 
the configuration. The minim-um drag of the wing was only a small part 
of the minimum drag of the complete configuration. 

6. The variation of drag with lift coefficient indicated that 

the resultant aerodynamic forces on the -wing and tail were approximately 
normal to the chord lines. 

7. The static stability of the complete configuration varied with 
lift coefficient at subsonic speeds. The somewhat irregular variation 
of static stability with Mach n-umber is apparently an effect of the 
wing on the configuration. 

8. The damping factor for the wingless configuration was practi- 
cally invariant with Mach n\miber. Addition of the wing caused the 
damping factor to have rather large variations with Mach number in the 
transonic region. 

9. The converging rear portion of the fuselage caused a downflow 
over the tail which became greater as the Mach mmiber increased. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I 


FUSELAGE NOSE AND TAIL ORDINATES 



X >j 








' ■ 




X 

(in.) 

(to.) 

0 

0.168 

0.060 

.182 

.122 

.210 

.2lv5 

.224 

.480 

.294 

.735 

.350 

1.225 

.462 

2.000 

.639 

2.450 

.735 

4.800 

1.245 

T.35O 

1.721 

8.000 

1.849 

9.800 

2.155 

12.250 

2.505 

13.125 

2.608 

14.375 

2.747 

14.700 

2.785 

17.150 

3.010 

19.600 

3.220 

22.050 

3.385 

24.500 

3.500 



TABLE II 


MODEL MASS CHARACTERISTICS 


Model 


Weight 

(lb) 

ly - 

(slug ft^) 

Center- of- gravity 
location 
(percent M.A.C . ) 

126.0 

8.91 

12.4 

107.0 

10.87 

-4.6 

88.8 

10.03 

12.3 



Side view 



Rear view 


Wing 

Aspect ratio 3.00 

Tat>er ratio 
Total area^sqff 2.68 
Exposed dreasg ft I .94 
Dihedral^ deg .00 

M.A.C.^in. 12.03 

Horizontal tail 

Aspect ratio 3.00 

Taper ratio 

Total drea^ sg ff .668 
Dihedral, deg -00 

M.A.C., in. 6.01 

Vertical tail 

Exposed area, sqft I I ' 


Figtire 1.- General arrangement of models A and B. (All dimensions are 

in inches.) 
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Vertical tail 

Exposed area 2.22 sq ft 


Horizontal tail 
Total area 
Aspect ratio 
Dihedral 
M.AC. 


.668 sq ft 

3.00 

.0 deg 

6.01 in. 


Figure 2.- General arrangement of model C. (All dimensions are in inches.) 
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L- 59511 



L- 59 512 



L-59515 


(a) Model A. 


Figiire 3*- Models tested. 
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(b) Model C. 



L- 61^-9 6 6 


Figiire 3-- Concluded 
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Figure k.- Tail surfaces^ model C. 



l-61j.6o3 



Figure 5*- Model for determining drag of vertical tail. (All dimensions 

are in inches.) 
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(a) Model A. 



L- 59 675 



(b) Model C. 

Figure 6.- Models on launcher. 



L-65550 
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Figure J.- Measured torsional rigidity of wings. 



Figure 8 .- Static pressvire ratio. 
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Figure 9»- Reynolds niimber of tests based on mean aerodynamic chord. 
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M=0.84 M=1.04 M=1.19 11=1.25 


a, deg 


(a) Model A. 
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□ O a decreasing 
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(b) Model B. 



M=1.23 M=1.34 

a, deg 
(c) Model C. 


Figure 10.- Variation of lift coefficient with angle of attack during 

several typical oscillations. 


4o 


MCA RM L50KL6 



.7 .8 .9 1.0 1.1 1.2 1.3 1.4 1.5 

M 

(a) Models A and B. 



(b) Model C. 


Figure 11.- Lift-curve slopes of complete models. 
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Figure 12.- Lift-curve slopes of wing and tail. 
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Figure 13 •- Lift effectiveness of elevator. 
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Figure l4.~ Summary of buffet and maximum lift information. 6 = -2.0°. 
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Figure 15 * - Minimum drag coefficients of complete models. 



Figure 16.- Minimum drag coefficients of model components. 
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(a) Model A. 
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(b) Model B. 


Figiire I7.- Effect of lift on drag 
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Figure 18.- Maximum lift-drag ratios. 
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Figure 19.- Lift coefficients at which maximum lift-drag ratios occ\ir. 
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(a) Model A. 
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(c) Model C. 


Figure 20.- Periods of oscillations. 
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Figiire 21.- Static stability of complete models. Center of gravity at 
12.4 percent of the mean aerodynamic chord. 



M 


Figure 22.- Static stability increments caused by wing and fuselage. 
Center of gravity at 12.4 percent of the mean aerodynamic chord. 


Aerodynamic center, percent M.A.C. 
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Figure 23.- Aerodynamic -center location 


Sap jad 


NACA RM L50KL6 


k9 


.8 


^1/2 


7 


















V c 
I>c 

/T 

) 

).04 





















— Model A 
■ — Model B 

— Model C 









Vv 

> 




— 













1 

! 














O Cl«0 
□ Cl»0.4 



t> 






o Cl«0.4 
A Cl~“*2 


" 








--A 























.7 .8 .9 1.0 1.1 1.2 1.3 1.4 1.5 

M 

(a) Time to damp to one-half amplitude. 
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(b) Pitching-moment damping factor 
Figure 24.- Damping characteristics of short-period oscillations. 
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(a) Trim lift coefficients. 
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(b) Trim angles of attack. 

Figure 25 .- Longitudinal trim characteristics. 
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Figure 26.- Effectiveness of the elevator in producing pitching moment. 



Fig\ire 27.- Change in trim lift coefficient per degree change in elevator 

deflection. 
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Figure 28.- Pitching -moment coefficient at zero angle of attack and zero 

elevator deflection. 
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Figure 29-- Downwash angle at zero angle of attack. 
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